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ry | \ The dynamics of all states of a qubit system with arbitrary, even time 
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dependent, one qubit Hamiltonians and two qubit interactions is realized as 



' one and the same diffusion process for systems with time dependent statistical 

OO ' 

weight in a space whose dimension grows only linearly with the number of 
qubits. The ensuing Fokker-Planck equation for the corresponding nonpos- 

o ' 

itive probability density is equivalent to the von Neumann equation for the 

Q_i| quantum state. Presumably the effectiveness of the stochastic process as a 

+-> ' 

numerical quantization method decreases as the number of qubits grows. 
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In the past few years, theoretical and experimental research on quantum computation, 



mainly sparked by the discovery of the quantum factorization algorithm by Shor [|J], led 
to a growing interest in quantum models with finite dimensional state spaces. This 
result confirmed a former foresight of Feynman j|] who had anticipated that amazing com- 
putational capabilities were inherent in quantum systems due to the exponential growth of 
their state space dimension as a function of their physical dimension. || Nevertheless, in 
spite of a general feeling, there is no formal proof that classical stochastic systems have 
a less effective scaling of their computational power versus physical dimensions . Apart 
from the relevance for the foundations of quantum mechanics, this issue makes the quest for 
classical stochastic models of finite dimensional quantum dynamics particularly appealing. 
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It was still Feynman who suggested that if one wants to get rid of quantum amplitudes 
and construct a quantum theory only in terms of probabilities, one has to use generalized 
non-positive probabilities. In fact a local traditional probabilistic theory always obeys Bell's 
inequalities []7|, contrary to what happens for quantum mechanical states. The aim of this 
paper is to prove that quantum evolution of qubit systems, namely arrays of two state 
quantum systems, with two qubit interactions, is equivalent to a suitable classical stochastic 
dynamics of particular non-positive probability densities in a phase space whose dimension 
is proportional to the number of qubits, which is the crucial feature of a classical state space. 

To fit the issue in the due context, it is worthwhile citing that Moyal |§ and Wigner || 
showed long ago how quantum states of point particles can be represented as nonpositive real 
functions on classical phase spaces. Their dynamics, which in the classical limit obeys the 
Liouville equation for the probability density, is nonlocal in space, i.e. their evolution cannot 
be described in terms of a partial differential equation of finite order. As a consequence its 
description as some sort of generalized classical diffusion process is not viable. 

If instead one considers, to begin with, a system of two interacting qubits, their generic 
(mixed) state can be expressed as an at most quadratic function of the Pauli generators 
of the corresponding sl(2,C) algebras. If these generators are replaced by real variables, 
these states can be interpreted as nonpositive probability densities giving rise to expecta- 
tion values coinciding with the quantum ones. On the other hand their quadratic character 
allows for the description of their time evolution by means of a second order partial dif- 
ferential equation. It is shown that this equation can be read as the master equation of 
a suitably generalized diffusion process. As to the generalization involved, if one wants to 
define statistical ensembles for non-positive probability densities, one has to introduce for 
each element in the ensemble a real statistical weight, which, as it is only natural to expect, 
will itself evolve in time according to a suitable differential equation. Although presented 
for notational simplicity with reference to a two qubit system, the result is general and the 
phase space is always a direct product of single qubit phase spaces. 

Stated in these terms the result seems to prove that quantum computational systems 
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can be effectively simulated by classical stochastic ones. Actually, since, as explained below, 
effectiveness seems to be achieved only for short computational times, the problem is subtler 
and deserves further investigation. A still deeper distinction between true quantum compu- 
tation and stochastic simulation in its present form resides in the lack within the latter of 
a mechanism of wave function collapse. This of course does not hold for NMR computing 
T0| . |TT1] where no wave function collapse occurs. 



Before proceeding to the construction of the quantization scheme it is worthwhile re- 



marking that it has nothing to do with Parisi-Wu quantization fl2|| , which is based on the 
addition of stochastic terms to a fictitious dissipative dynamics which has the ground state 
as a fixed point. On the other hand, while Nelson's approach ]13[ deals with real time evo- 



lution, it is tailored for point particle dynamics and it does not follow Feynman's suggestion 
about nonpositive probabilities and Moyal-Wigner viewpoint about treating pure and mixed 
states on the same footing. 

To be specific, consider a generic two qubit system formally identified with a two spin 
1/2 system, whose Hamiltonian, apart from an additive constant, has the traceless form 

H = -B^S? - S]Jij Sj, (1) 

where summation over repeated indices and apices is implied, S},Sf denote the two spin 
1/2 vector components defining the usual Clifford algebra 



°i > °j 



^a^ijkS^, S^Sj — -5 a/ 35ij, (2) 



while B 1 , B 2 and J denote possibly time dependent c- number vectors and a 3 x 3 matrix. 
The most general state, i.e. unit trace non negative Hermitian operator, reads 

p = ± (l + 4 (s?) St + 16 (SIS') SIS') , (3) 

where 1 is the unit operator in the four dimensional Hilbert space. Dynamical evolution is 
completely characterized in terms of the expectation values appearing in Eq. (H), i.e. by 
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' = -i ([ S i> H \) = ~ £ ijk (B) (pi) - Jji (SlS? 
d(S?) ' r SlH]) = -e m (B) (Si) - J, (S}St)) (4) 



dt 



J I Cl C2 

' J'/ _ _ r1/q1q2\ _ r 2/q1q2\ £imfc t / M\ £jmk j I A 2 



where in the last equation anticommutation properties of the Clifford algebra (0) were used. 

If in Eq. (|j) spin components outside expectation values are promoted to the components 
S}, Sf of ordinary real vectors S 1 , S 2 , then for the generic traceless observable one gets 

/ ds 1 j ds 2 (^l + 9 -^4^-)p(s\s' 
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a" a "lc2\ ^P 1 !^! r 2 <\S 1 \<R 2 

Cj Sj + CijS; Sj ) = - — y-=i ~ r (5) 

J/ J dS 1 J dS 2 p(s\S 2 ) 

r-i < I 1 r 2 <\s 1 \<R 2 

for arbitrary values of r 1? i?i,r 2 , R 2 , by which p can be interpreted as an unnormalizable 
non-positive probability density evolving according to the master equation 

do d(S?) d(S}S 2 ) 

op = A^ls? + 4 V ; 3 ' SIS 2 . (6) 
dt dt dt 1 3 w 

Here and henceforth only polynomial probability densities obtained from Eq. (|3]) are con- 
sidered and called quantum states. Using Eq.s (II), Eq. (^) is seen to be equivalent to 

°P fZa a ffa\ 9 P (ol, 9o2\ d P (o2 A jTol\ °P 



= (B a A S a ) ■ —S- (S l A JS 2 ) ■ —4 (S 2 A J T S . , 

dt > dS a v ' dS 1 v 7 dS 2 

l-i d -dp 1-n d it dp 
+ -S 1 A — — ■ J— b- + -S 2 A — — ■ J T — b- (7) 
4 as 1 <9S 2 4 OS 2 dS 1 

° Cl 9 P , fc?l a 7C2\ 5 c2 



v y ^ 2 v > dS 1 dS 2 ' 

where standard vector notation was used. While the first three terms above give the Liou- 
ville equation for two interacting classical tops, the following ones cannot be interpreted as 
resulting from a deterministic dynamics in the space of S 1 and S 2 since they contain second 
partial derivatives. In the following it is shown that the above equation can be interpreted 
as the master equation of a generalized diffusion process in this space. 

The generalization implied consists in considering a statistical ensemble of equivalent 
systems each one endowed with a nonpositive real statistical weight w obeying an ODE 



dw = h(z, t)wdt (8) 
and whose states are represented by coordinate vectors z obeying a Langevin equation 

dz = f(z, t)dt + 9j(z, t)dujj, (9) 
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where the gaussian random variables duji have zero mean duji = and covariances defined 



by u>i(t)d)j(t') = 5ij5(t — t') : or less naively by 

rt+dt _ rt'+dt' 

[Ui(t + dt) - Ui(t)] [uj(t' + dt') - Uj(t')] = J dt J t dt' S(i - t'). 

A statistical ensemble (z, tUi)i=i,2,..iv defines a nonpositive probability density p(z, t) 
whose averages coincide with the weighted ensemble averages 

1 N r 1 N 

P(z,t) = —^Wi5{z-Zi(t)) =>- / dzpk(z) = — ^Wik(z,i) (10) 

i=l i=l 

and whose time evolution contains one term more than usual: 

I N q l N Q 

dp = ~J^Yl W i d ^ ' ~Q=$> (Z- Zi(t)) + JjYl 5 (Z~ dnJ i = ' ( d *P) + hpdt. (11) 

Starting from this equation one can obtain the generalized Fokker-Planck (F-P) equation: 

~k = ~d~z' t)p ) + 2^d~z®d~z' ^ t] ® t)f>] + H ^ )P - (12) 

Of course one could even adopt a more traditional, though less insightful viewpoint, by 
treating w as a further dynamical variable and introducing a positive probability density 

1 N r+oo 

p P (z,w,t) = — J25(z- Zi(t))5(w-Wi(t)) p(z,t) = / wp P (z,w,t)dw, (13) 

iV i=1 J-oo 

in which case the additional term in the generalized F-P equation is simply the transcription 
of a term in the ordinary F-P equation for pp(z,w,t): 

/+oo Q r+oo 
w——[h(z,t)wpp(z,w,t)]dw— / h(z,t)wpp(z,w,t)dw — h(z,t)p(z,t), (14) 
-oo OW J-oo 

where it was assumed that w 2 p P (z,w,t) — > as w — > ±oo. 



While this alternative setting leads to the usual form of the F-P equation, the price to pay 
consists, when evaluating averages, in replacing the generic dynamical variable k(z) with the 
product wk(z), which is of course the remnant of the nonpositive character of the original 
probability density and then the source of the possible violation of Bell's inequalities. 

To be specific, if one defines the vector fields 



F 1 = IS 1 A J' 
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(15) 



one can show that Eq.(|7|) is the F-P equation of the system of stochastic differential equations 



dS 1 



-B 1 A S 1 + 5S 1 A JS 2 + F 1 + 6 j^j 2 S 1 - -S 1 
+ ^ A S 1 + j^ + d^S 1 AJ^ + d^ s 1 s 1 



dt 



dS 2 



2 2 ' - ' s 2 

-B 2 A S 2 + 5S 2 A JTS 1 + F 2 + 6 
dff 



(16) 
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2 



s 2 



s 2 -^s 2 



dt 



S 1 
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dw = w 



' f ° .^ + -4-.^ +15 



2 XdS 1 



OS 2 



S 1 



+ 



s 2 

2\ 3 



dt, 



(17) 
(18) 



where ff , r\ , £l,£2 are zero mean independent gaussian ^-correlated random variables: 



(19) 



The contribution to the F-P equation due to terms in ff 1 , apart from operators vanishing on 
quantum states and before averaging on ff 1 , reads 



1 

dt 



' 9 r- ■ (dff 1 A S 1 ) • (dff 1 A S 1 ) + \^=r ■ (dff 1 A S 1 ) dff 1 ■ J- 9 



SdS 1 v J dS 1 v ' ^dS 1 v ' 7 ' OS 2 

where here and henceforth the symbol <— d — > denotes that the corresponding derivative 
operator applies even to the left, and after averaging it reads 

4#.4*-i.* + i(*A4,y-?-t. 

4 OS 1 4 dS 1 4^ dS 1 J OS 2 

Of course analogous relations hold true for terms in ff. 
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As to the terms in and £ 2 , after averaging they give rise respectively to 
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for I; 1 and an analogous expression for £ 2 . 



S 2 + 3 



s 2 - 



d_ 

dl 2 



Adding these terms to the Liouville operator of the deterministic terms in Eq.s (|16|,lT) 
and to the ones generated by Eq. (p^8|) for the statistical weight according to Eq. fll2|), one 
gets the proof of the equivalence between the stochastic differential equations (p~6| , p~T| . p~8|) and 
the F-P Eq. (g). 

Of course the present result and its proof hold basically unchanged for systems with more 
than two qubits, provided that there are only pair interactions. To be specific, consider, in 
its proof, what turns out when taking account of the time evolution due to a particular two 
qubit interaction. The crucial point is that, if p is an arbitrary polynomial at most linear in 
every spin vector, the J dependent part of the right hand side of Eq. (^) coincides with what 
is got considering spin components as operators and taking the corresponding commutator 
appearing in the von Neumann equation. Then in general to obtain the master equation one 
simply has to introduce operators of the same form as in Eq. ([7]) for each pair interaction, 
which obviously leads to an immediate generalization of equations (|16|,[17],|18|) . 

Some comments are in order to the numerical implementation of the quantization 
method. Once fixed an initial quantum state, one takes an ensemble sampling the corre- 
sponding probability density for S a < M where M is properly chosen. Than the numerical 



simulation will work for 



< m < M until the effect of the lacking sampling for 



S c 



> M 



will be felt at 



< m. Of course the reliable simulation time will be the longer the greater 



M and the lower m are chosen. In other words finite statistical ensembles undergo an aging 
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effect. So that if one wants to simulate for longer times, one has to reset them periodically. 
In doing this at each step a fraction of the statistical information encoded is lost. 

A thorough numerical analysis is needed to assess the real power of the method as 
the number of qubits grows. A way out of the seeming contradiction between the linear 
growth of the simulation space in the present quantization method and the exponential 
growth, with the number of qubits, of the computational resources needed for simulation 
of quantum systems @] in the author opinion is to be found precisely in this aging effect. 
An indication in this direction comes from the observation that the logarithmic derivative 
of the statistical weight is proportional to the number of interaction terms, which is to be 
interpreted as a growing compensating effect for the depopulation and overpopulation effect 
of time evolution in different regions of the S a space as the number of pair interactions grows. 
This would produce an aging effect even for an hypothetical infinite statistical ensemble 
sampling the whole space, as can be guessed for instance for the totally disordered state, 
whose density matrix is proportional to the identity matrix. While the probability density 
is obviously constant for this state, this is not the case for the statistical ensemble. In 
particular the ultimate fate of every stochastic trajectory is to run away to infinity, due to 
the cubic deterministic term. While this term can be eliminated with alternative forms of 
the stochastic differential equations, probably the present one is the most convenient as to 
the behavior close to the origin. 

It is relevant to remark that the aging of the statistical ensembles reconciles another 
apparent contradiction inherent in the use of a typically irreversible diffusion process to 
describe reversible unitary quantum evolution. 

It should also be stressed that, although spin notation in this paper played in general 
only a formal role, the presented result applies to spin systems like arbitrarily disordered 
anisotropic Heisenberg systems. In the case of long range interactions of course the quan- 
tization method is numerically less effective due to the quadratic growth of the number of 
stochastic terms. 

Finally it should be remarked that the present setting has a natural generalization to 
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quantum open systems in the Markoff approximation in terms of additional stochastic terms. 
One can naturally introduce a fluctuating part in the B a vectors, which simply amounts to 
the addition of further terms in the F-P equation. 
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